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Vector modulational instability induced by parametric resonance in periodically 
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We study the modulational instability induced by periodic variations of group- velocity dispersion 
and nonlinear coefficients in a highly birefringent fiber. We observe, for each resonance order, the 
presence of two pairs of genuine vector type sidebands, which are spectrally unbalanced between 
the polarization components for nonzero group-index mismatch, and one pair of balanced sidebands 
emerging and dominating at increasing group-index mismatch. As the conventional modulational 
instability manifests itself, it is partially suppressed by the proximity of these new unstable regions. 

PACS numbers: 42.65.Sf, 42.65.Ky, 42.65.Tg 



I. INTRODUCTION 

In classical mechanics parametric resonance (PR) is a 
well-known instability phenomenon which occurs in sys- 
tems the parameters of which are varied periodically dur- 
ing evolution [IJ[2]. For example, a harmonic oscillator 
the frequency of which is forced to vary in time will be- 
come unstable if its internal parameters and the ampli- 
tude of the frequency variation happen to be inside spe- 
cial regions, known as resonance tongues. The study of 
the properties of resonance tongues has a long history 
and relies on a variety of geometrical approaches [3] 0] . 

It is natural that such a general phenomenon was as- 
sociated to the equally important instability process that 
is ubiquitous in infinite dimensional dynamical systems: 
modulation instability (MI), also known as Benjamin- 
Feir instability [5]. MI is known to exist in different 
branches of physics such as fluid-dynamics [6], plasma 
physics [7], Bose-Einstein condensates [8] and solid-state 
physics |5]. In nonlinear optics [TO] , it manifests itself 
as pairs of sidebands exponentially growing on top of a 
plane wave initial condition, by virtue of the interplay 
between the cubic Kerr nonlinearity and the group ve- 
locity dispersion (GVD). In optical fibers it leads to the 
breakup of a plane wave into a train of normal modes of 
the system, i.e. solitons [TT] . 

The link between PR and MI has been established in 
relation to the periodic re-amplification of signals in long- 
haul telecommunication optical fiber cables [12]. This 
was based on a nonlinear Schrodinger equation (NLS) 
where the coefficient of the nonlinear term is varied 
along the propagation direction. Importantly, this pe- 
culiar type of MI occurs in both normal and anomalous 
GVD. This prediction was later partially verified in ex- 
periments, see [13] . 

Moreover, in long-haul fibers, dispersion management 
is a commonly used technique which introduces periodic 
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modulation of fiber characteristics. The possibility of in- 
stability phenomena disrupting adjacent communication 
channels has been thoroughly analyzed, see e.g. [T4HT6] . 
Specifically, in [14] the partial suppression of the con- 
ventional MI in anomalous GVD due to a large swing 
dispersion management is discussed, while in [15] the de- 
generate case of zero average dispersion was studied. The 
combination of both loss and dispersion compensation is 
studied in [16] . The main interest in those works was on 
step-like variations of the GVD coefficient. 

At the same time the effects of smooth periodic or 
random variations of fiber parameters were studied in [T71 
H8] . Also some work has been done on the effect of the 
perturbation of fiber parameters on soliton propagation 

It turns out that the variation of dispersion and non- 
linearity can enhance or suppress the PR, while higher 
order nonlinear effects such as self-steepening proves less 
important. Quite surprisingly, experiments on micro- 
structured fibers have been reported only recently for the 
first time, see Ref. [20], where a photonic-crystal fiber 
(PCF, [21]) of varying diameter is used. In that experi- 
ment, the dispersion is periodically switched from normal 
to anomalous, but this feature is not required to achieve 
PR, while the effect of Raman scattering plays an impor- 
tant role in the relative magnitude of the PR peaks. 

The conventional explanation is in term of a grating- 
assisted phase matching process [12] [16j [20] , but it was 
verified in Ref. [22] that this approximation is inaccurate 
if the period of parameter variation is comparable with 
the length scale at which the nonlinear processes occur. 
In Ref. [22] it was proved that an accurate description 
must be based on the Floquet theory [2, 3 and the use 
of regular perturbation techniques, such as the method 
of averaging [23 . 

The study of birefringent fibers permits to observe a 
variety of new physical phenomena, which are ascribed 
to the presence of cross-phase modulation (XPM) terms 
[24] . The MI in birefringent fibers (vector MI) occurs also 
for normal GVD and was extensively studied in the past 
not only in highly birefringent fibers (HBFs) [25-27 but 



also for weak birefringence [28] . The effect of a step- wise 
variation of birefringence was considered in Ref. [29] , in 
the weakly birefringent regime, while highly-birefringent 
fibers with step- wise variations of dispersion were stud- 
ied in Ref. [18], in a dispersion- management scenario of 
alternating GVD sign. This two last works apply rig- 
orously the Floquet theory, but ignore completely the 
group-index mismatch. 

The possibility of tailoring the PCF birefringence, see 
[2T] and references therein, and of obtaining a smooth 
reproducible profile of fiber parameters by advanced fab- 
rication techniques, [20] , permits to achieve PR instabil- 
ities on a short distance and to explore different birefrin- 
gence regimes and the effect of group velocity mismatch. 

In the present work we study parametric instabilities in 
a HBF with varying GVD and nonlinear coefficient. We 
provide accurate analytical estimate of PR peak detuning 
and gain and contrast them to the numerical application 
of Floquet theory and to split-step simulations. We ob- 
serve the existence of two families of MI peaks at each 
PR order: one exhibits a behavior similar to conventional 
(i.e. with constant parameters) vector MI, while the other 
resembles scalar MI and is the dominant MI process for 
large group-index mismatches. Finally we found that, at 
large group-index mismatch, the conventional vector MI 
is partially suppressed for large variations of parameters. 



II. MODEL EQUATIONS AND ANALYTICAL 
ESTIMATES 

A. Incoherently coupled NLS and linearized 
equations 

According to the conventional approach, [24 , the prop- 
agation in HBF can be described by two incoherently- 
coupled NLS equations (ICNLS), which read as 



id z A, ±i-d t Aj 



1 



; z ^j 



/3 2 (z)d tt A j + 



7(*)(l^ 



B\A 3 .j\ 2 )Aj = 0, with j = 1,2, (1) 



where /3 2 and 7 are normalized GVD and nonlinear co- 
efficients, p 2 (z) = P 2 (z)/73 2 and 70) = 70)/7 ; £2(2) 
and j(z) are the physical GVD and nonlinear coeffi- 
cients, respectively, and the superscript denotes their 
mean values. 7 and ($2 are assumed to be equal for 
the two polarizations and periodic functions of z. Fi- 
nally z = Z/Z n \ is the dimensionless distance in units 
of the nonlinear length Z n \ = (7°P^) _1 , and t = (T — 

(vg) Z)/T s is the dimensionless retarded time in units 



of T s = y Z n \\f3 2 \, v^ is the mean group velocity, and 

8 — Z n \/T s [(^ 1 )i — (Vg 1 )2\ is the normalized group- 
index mismatch between the two polarizations. P t is the 
total input power injected in the fiber, and A\$ are the 
dimensionless slowly varying modal amplitudes of the two 



polarization components scaled by \[Pf The XPM coef- 
ficient B is used throughout the paper since the ICNLS 
model can be applied to other physical settings [30 . In a 
HBF, we can neglect the phase mismatch term due to co- 
herent coupling and set B = 2/3, thus we have that Ai^ 
correspond to the mode polarized along the fast and slow 
axis, respectively. 

We look for for a steady state solution of (HI) in the 
form Ai 5 2 = \fP1~2 exp {i(j)i^{z)): it can be verified 
that </> li2 (z) = (Pi,2 + BP 2 ' 1 )f^ oo j(z f )dz f . We then 
perturb this steady state by adding a small complex 
time dependent contribution a\^{z, t), i.e. Ai i2 (z,t) = 
(y/Pi,2 + £Q>i,2(z,t)) exp M>i, 2 (z)), with £«1. Insert- 
ing this Ansatz in Eq. (JTl) and taking only the terms 
which are first order in £, one finds that a\^ obeys the 
following equation: 



id z CLj ± i-dtdj - -p 2 (z)d^aj + 7(2) [Pj(aj + a*)- 



By/Pj^~(a 3 . 



j+ a 3-jJ 



0, j = 1, 2. (2) 



We further assume that the input light is polarized at 
an angle of 7r/4 with respect to the fast axis, i.e. P = 
Pi = P2 = 1/2, which significantly simplifies our calcula- 
tions and that GVD and nonlinearity exhibit the simplest 
possible periodic behavior 



A>0) = A) + P(z) = A) + hpi cos Az, 
7(2) =70 + 7(2) =7o + fryi cos Az, 



(3) 



where generally /3q = +1 for normal (anomalous) GVD 
and 70 = 1; A is the normalized spatial angular frequency 
for the parameter oscillations. The forcing amplitude 
is controlled by the parameter h, which must be small 
to guarantee the validity of our perturbative expansions. 
However, we find below that our estimates are reliable 
even for h ~ 0.5. Finally we substitute in (|2| the Ansatz 

aj (z,t) = af(z)e- iut + af(z)e iujt , j = 1,2, 

which permits to cast the linearized system in the form 
of a 4 th -order linear ODE system 



dz l 



where 



H(z)\(/>), H(z) 
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(4) 



\(j)) = (ui,Vi,U2,V2) T , 

Uj = a? + a/*, Vj = dj - af*, j = 1,2 (5) 

and we defined c\{z) = — ^-/^(z) = c\ + ZiCiCosAz, 
c 2 {z) — ci(z) — 2j(z)P = C2 + hc2CosAz, b(z) = 
—j(z)BP = bo + hb cos Az, and v = isq = — |cj. By 
replacing pi in these definitions we can naturally split 



the Hamiltonian matrix H(z) into average and oscillating 
parts, i.e. H(z) = Ho + hH(z). 

Eq. Q can be rewritten as a system of two coupled 
Hill's equations, i.e. linear oscillators with periodic vari- 
ation of natural frequencies, but it is more practical to 
deal with the original first order system directly. 



B. Calculating position and gain of PR peaks by 
the averaging method 

We first present the relation which provides the values 
of PR detuning. We discussed extensively in [22] how to 
apply the classical theory of parametric resonance [2] to 
problems of instability in fiber optics involving varying 
parameters. 

Parametric resonance is a phenomenon which is ac- 
curately described by a relation between the natural fre- 
quency of the unperturbed oscillator and the forcing term 
frequency. Thus we have to impose that Hq has real 
eigenvalues, which in turn implies PR is incompatible 
with conventional MI. From a physical point of view this 
is justified by the fact that conventional MI is generally 
a much stronger instability effect. 

Since the matrix Ho is traceless [31] and given the 
choice of identical GVD and nonlinear coefficients and of 
the particular polarization state, we can write its eigen- 
values as ±Ai and ±A2, with 



Al2 — 



c l c 2 
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A single parametric oscillator is destabilized if the un- 
perturbed system oscillates at half an integer multiple of 
the forcing frequency. In the present case we have two 
coupled oscillators and the scenario is more complicated. 
We must consider four independent conditions: 

2Ai, 2 = raA, (7) 

which we denote as vector MI band (V-band) and 

Ai ± A 2 = raA, (8) 

denoted by scalar-like MI band (S-band), where m is the 
PR order and the reason of the definitions will be made 
clear below. In each case we obtain a polynomial in uj^ 
the zeros of which gives us the PR detuning. The two 
polynomials are reported in Appendix [Aj 

The relations between the spatial frequency of external 
forcing and the eigenvalues of Ho can also be obtained 
by the method of averaging [23] . In its simplest formula- 
tion it is based on the method of variation of constants 
for inhomogeneous differential equations. This in turn is 
equivalent to transforming the system of Eq. Q to the 
interaction picture, i.e. the evolution of the slow variables 
is governed by 



= e iH z\ 



The averaging process is used to eliminate the remain- 
ing oscillating terms from Eq. Q . In the right-hand side 
we find elements with spatial periods obtained by linear 
combinations of A and Ai } 2, which are in general incom- 
mensurable. Thus the method of averaging needs to be 
generalized by performing the integration over an infinite 
range, i.e. 

id z \4>) I = h(H I (z))\4>) I , 

(Hj(z)) = lim \ f H I (z f )dz f (10) 

Z^oo Z Jo 

It is clear from this matrix expression how to obtain 
the four PR conditions, since the above-mentioned res- 
onances correspond to the presence of nonzero average 
elements in the interaction Hamiltonian. Most impor- 
tantly we can estimate the peak gain, at first order in h, 
by solving for the complex eigenvalues of the averaged 
interaction Hamiltonian. 

As in the scalar case of Ref. [22] , the first order aver- 
aging method provides us with an estimate of the peak 
gain of the l st -order PR, which are reported in Appendix 
|B] In order to estimate the gain of higher-order PR, a 
higher-order perturbation theory is demanded, but this 
is outside the scope of this work. 

In the next paragraph we present the numerical char- 
acterization of the PR phenomenon, in the from of res- 
onance tongues and output spectra of split-step simula- 
tions and compare it to our analytical estimates. 

III. RESULTS AND DISCUSSION 

Throughout this paragraph we set fio = +1, normal 
GVD, A = 10 and 71 = — fi\ = 1, the latter associ- 
ated to the maximum gain in the scalar case [22]. As 
a guide for our considerations we study first how the 
properties of PR sidebands as a function of the group- 
index mismatch S. This parameter was neglected in 
the past [18 , but it plays here a crucial role. For 
our choice of parameters, the conventional MI occurs at 
S > ^4/^070^/3 = ^2jl « 0.82. In Fig. fl] we report 
the analytical estimates as a function of S 01(a) the PR 
and ( |A2| )], (b) their respective gain 
)] (c) and the Stokes-antiStokes im- 



detuning [Eqs. (LAI 
[Eq. ([Bl]) and |B2 
balance of sidebands, defined as R\ 



dz ] 



) I = hH I (z)\$) I , withtf 7 



5 iHoZ H(z)e- iHoZ . 



(9) 



l a i / a i I? which 
is obtained by the eigenvectors of the averaged Hamilto- 
nian. The imbalance is defined only for one polarization 
mode, as for the other polarization component is exactly 
the inverse, on account of the conservation of total mo- 
mentum of Eq. . 

In Fig. [T] we observe that the S-bands occur at con- 
stant detuning between a pair of V-bands (a) , they have 
finite gain only if S 7^ (b) and are spectrally symmet- 
ric around the pump frequency [see (c)], which justify our 
definition of scalar-like bands. At around 5 ~ 3.5 the con- 
ventional MI unstable sideband crosses the S-band and 
the latter switches from plus to minus sign in Eq. ([8]), 
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FIG. 1. (Color online) Characterization of 1 st order PR as 
a function of normalized group- index mismatch (8), for pa- 
rameters varying with A = 10. (a) Resonant detuning; (b) 
Gain slope (gi/h) values (c) Unbalance of V-bands around uj 
for the fast axis; the slow axis exhibits the opposite behav- 
ior. The S-bands are always balanced and are reported in the 
region of nonzero gain. The following line convention applies 
to every panel: blue solid line is the V-band corresponding to 
A = 2A2, while the red dashed to the other V-band, A = 2Ai. 
The green dash-dotted line denotes the S-band A = A2 + Ai; 
finally the black dotted line corresponds to the other S-band, 
A = A2 — Ai. In panel (a), we also report the instability 
margins of the conventional XPM MI in the absence of per- 
turbations, as pink dashed line marked by asterisks. 



as can be noticed by carefully observing the range where 
the gain is zero in (b). Moreover their gain is constant 
over a wide range of 5. 

The V-bands amplitudes are perfectly symmetric 
around uj = for S = while they develop an asymmetry 
for 5 ^ 0: thus this PR bands have the same character of 
the conventional vector MI bands in the ICNLS system 
[25j |30j [32] and this explains our definition. Finally they 
are increasingly split apart as S increases. 

The brightest, i.e. largest gain, peak is for 6 < 0.52 a 
V-band, then, for 0.52 < 6 < 2.5 an S-band. We will 
discuss below what happens beyond 5 « 2.5, where a V- 
band exhibits a gain larger then the S-band: the numer- 
ically computed resonance tongues shows a complicated 
structure where conventional MI and high-order PR co- 
exist and the gain predictions prove inaccurate. 

We then discuss the structure of instability tongues for 
4 different cases: (i) 5 = in Fig. [2J (ii) 5 = 0.4 in Fig. 
[3] (iii) S = 1.15 in Fig. g] and (iv) S = 4 in Fig. |5J 

In Fig. [2] we observe two V-bands, the first with a large 
peak gain and the second with a much weaker gain. De- 
spite we consider S = as in [18] , we have the important 
difference that, here, the GVD varies smoothly and is al- 
ways in the normal region, instead of the step-wise with 
alternating sign presented in that paper, so that we do 



FIG. 2. (Color online) Resonance tongues for 1 st order PR, 
with A = 10 and 6 = 0. The color scale corresponds to the in- 
stability gain The dashed lines denote the predicted positions 
of PR peaks, while the corresponding cross-marked lines rep- 
resent the numerically obtained position. The corresponding 
maximum gain is showed in the inset (solid lines) as a func- 
tion of the perturbation strength h and is compared with the 
analytical predictions (dashed lines). The line colors in the 
inset correspond to those used in the uj — h diagram for the 
peak detuning positions. 




FIG. 3. (Color online) Same as Fig. [2] with 5 = 0.4. 



not observe conventional ML Finally we report the posi- 
tion of the S-band, which exhibits vanishing gain. It can 
be verified numerically that the V-bands grow spectrally 
symmetric in both polarizations, as expected for 5 = 0. 

Figures |3] and [4] show two similar situations: the main 
difference is that in the former the lower detuned V-band 
are the brightest MI peak while in the latter the central 
S-band has overcome V-bands as the brightest gain side- 
band. Conventional MI occurs at small uj for S = 1.15, 
but is not reported since it is not influenced by PR and 
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FIG. 4. (Color online) Same as Fig. [5] with S = 1.15. 
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FIG. 5. (Color online) Same as Fig. [2] with S = 4. The red 
line shows the gain curve of conventional vector MI, which is 
partially suppressed at large h. 



is almost independent of h. Before concentrating more 
on the case of Fig. [3J we finally present, in Fig. [5j the 
resonance tongues at large detuning S = 4. We observe 
that the sideband structure of PR is still dominated by 
the PR S-band, while higher order sidebands [at uj « 2.8 
(V), uj w 4.4 (S) of 2 nd order and at uj w 5.4 (S) of 3 rd or- 
der] are interleaved with the 1 st order ones. Moreover the 
conventional MI and the V-band at uj = 5.6 are partially 
suppressed for large ft, on account of the proximity of the 
higher-order peaks. At such values the first order esti- 
mate of gain is clearly inadequate (see inset), as in general 
occurs for the V-bands which coexist at large detuning 
with the conventional MI and higher-order PR peaks. We 
thus observe that for S > 3.5, where the conventional MI 
occurs beyond the brightest PR peaks of scalar-like type, 
the variations of parameters enhances spectrally symmet- 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




i t~*~t 1 * / * — 

Iff M f 

Ix j Ik! k 




* * r 

k f * 


J^r 


t t * 0.5 
xxx h 
i i k 



0.25 



10.15 



0.05 



4.2 



4.4 



4.6 



4.8 



FIG. 6. (Color online) 2 nd -order PR instability regions, same 
parameters and conventions as in Fig. [3] In the inset only the 
numerical results are reported, since analytical estimates are 
not considered here. 



ric, scalar-like, sidebands and suppresses the asymmetric 
sidebands which are commonly considered as the charac- 
terizing feature of MI in HBFs. 

This behavior is consistent with the suppression of vec- 
tor MI sidebands due to fluctuations of the fiber param- 
eters, which has for long precluded their observation in 
PCFs, see [33]. We tested the effect of periodic varia- 
tions of S and observe a reduction of the peak gain of 
V-bands for large average 5. Expressions of gain can be 
obtained, but are more involved than those presented in 
Appendix [Bj the PR detuning values are robust to this 
perturbations and the instability growth happens on a 
length scale larger than the period of the parameter vari- 
ations, so we decided not to explicitly consider variations 
of 5 here. 

We now complete the characterization in the case of 
Fig. [3j by showing in Fig. [6] the second order PR insta- 
bility regions for S = 0.4, which share the same features 
of the 1 st order ones, except the central S-band gain is 
already slightly larger than the smaller detuned V-band. 

Finally we include the output spectra obtained by solv- 
ing the system of Eq. (HI) by means of the split-step 
method, see Fig. [7[ We set all the parameters as above, 
ft = 0.9 and 6 = 0.4. 

We clearly identify the 1 st and 2 nd order PR, each of 
which is composed by three peaks. It is thus clear that 
the S-bands occur between a couple of V-bands and are 
symmetrically growing about the pump in both axis. The 
two V-bands exhibit an imbalance about uj = which is 
reversed from fast to slow axis on account of the conserva- 
tion of total momentum for the model of Eq. (IT]), compare 
panels (a) and (b). Moreover the two V-bands exhibit op- 
posite symmetry; consider the fast axis Fig.[7^a): the first 
peak at each order is characterized by the Stokes side- 
band outgrowing the anti-Stokes, while the second ex- 




tion at spatial (angular) frequency A, as in the scalar PR 
[22] : this is quite effectively explained by the theory of 
averaging. However there is a remarkable difference be- 
tween the two peaks: (a) the V-band involves only one 
eigenvalue (A2 = A/2) of H and grows upon a simple 
oscillation, while (b) the S-band involves both eigenval- 
ues of Hq and exhibits a beat of the fast oscillations, 
corresponding to Ai + A2 = A, and slow oscillations, cor- 
responding to A2 — Ai = A. 

Finally the imbalance as computed numerically from 
Fig. |8Fa) is smaller than that in Fig. flic); this a general 
trend: we observed that for large h the imbalance of V- 
bands is smaller than expected by the eigenvectors of the 
averaged interaction Hamiltonian. 



FIG. 7. (Color online) Output spectrum at normalized dis- 
tance z = 30 obtained by split-step numerical simulation. 
h — 0.9, 5 — 0.4 and the other parameters are as in the previ- 
ous figures; (a) axis 1 and (b) axis 2. We identify and classify 
the 1 st and 2 nd order PR peaks, as indicated by the text in 
the panel, by their spectral imbalance with respect to uj — 0. 
The small peaks near the pump components are the four- wave 
mixing product of the brightest l st -order V-band and S-band. 




IV. CONCLUSIONS 

In this paper we studied the effect of the periodic vari- 
ation of group- velocity dispersion and nonlinearity on the 
propagation of light in a highly-birefringent optical fiber. 
We showed that MI sidebands are effectively described in 
terms of parametric instabilities of a system of coupled 
oscillators and provide accurate analytical estimates of 
their detuning and gain. We considered only the nor- 
mal GVD regime and discovered the existence of two 
different kinds of unstable sidebands: the first, similar 
to conventional MI, appears as two pairs of sidebands 
which generally exhibit spectral imbalance around the 
pumps, while the second manifests itself only for non- 
zero group-index mismatch as a pair of spectrally bal- 
anced peaks: moreover the latter becomes the brightest 
unstable peak for large enough mismatch values. The 
vector sidebands and the conventional vector (XPM) MI 
are partially suppressed by the proximity of scalar-like 
PR sidebands. The PR peak position is widely tunable 
by varying the period of variation of parameters and the 
input power; this phenomenon could thus find interesting 
applications in quantum optics. 



FIG. 8. (Color online) Simulated evolution, on the fast axis, 
of the two main peak frequencies of 1 st order PR: (a) V-band 
at uj — 2.79 and (b) the S-band uj — 3.01. Red (blue with 
crosses) solid lines corresponds to the numerical evolution of 
(anti-)Stokes-sideband; dashed line are the average growth 
predicted by analytical calculations. 
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hibits the opposite behavior — the anti- Stokes dominates 
over the Stokes. This is analogous to what occurs in 
the proximity of zero dispersion due to the presence of 
higher-order dispersion [32] . 

The growth trend is presented in Fig. [8j for the two 
brightest first order peaks. The exponential growth of 
the unstable frequencies is superimposed to an oscilla- 



Appendix A: Estimate of resonant detuning: 
expression of polynomials 



In order to obtain the resonant detuning we substitute 
all the quantities defined after Eq. ^ in Eq. ([6]), then 
recast the conditions for 771-th order PR as a polynomial 



in cj^, which read as 



main points and results. Let us diagonalize Hq as 



P^ 2 ) = ^ 8 + ^(8^o 3 7o - 2/3 ^ 2 )c 6 + 
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for V-bands, Eq. (17]), and 
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id 



[S 2 A 2 (-m 2 ) - 4/3 o7 oPm 2 A 2 ] u 2 + ^-, 



for S-bands, Eq. ([8|. 



(A2) 



#o = VAV' 1 



where V is the matrix, the columns of which are the 
eigenvectors of Hq (not necessarily normalized) and A = 
diag[-A 2 ,-Ai,Ai,A 2 ]. 

The interaction Hamiltonian in Eq. ^ can be ex- 
panded as 



H T = VexpiiAz^^HVexpi-iAz^- 1 



but, since V does not depend on z, we can resort to the 
similar matrix 



H'j = expiiAz^^HV exp(-iAz), 



Appendix B: Method of averaging: instability gain 

The explicit calculation leading expressions of the PR 
peak gain is quite tedious, thus we summarize here the 



thus simplifying the resulting averaged matrix. 

The expressions of gain are different for each of the 
cases in Eq. ffl and |8]); we express them in compact 
form as 



aX 



4c?6 A 2 (A 2 - A 2 ) 



{[■ 



Aici + Pic 2 + C x b 



A 2 c x + B 2 c 2 + C 2 b 



y 



(Bl) 



for Eq. 0, with 2A 2 = A, with A x = b [±c\c%v(y + A 2 ) + (A| - A?)(A 2 + v) 2 } , A 2 

b [±c\c\v(v - A 2 ) + (A 2 - A 2 )(A 2 - v) 2 ] , B 1 = b c\ [(A 2 + 2v) 2 - A 2 ] , B 2 = b 4 [(A 2 - 2vf - A 2 ] , C x = 46 2 (c°) 3 - 

44 - (A 2 + vf\ (44 - A 2 + iv 2 + 2\ 2 v), C 2 = Ab 2 (c°) 3 - c? [44 - (A 2 - vf\ (44 - A 2 + Zv 2 - 2\ 2 v). 

I 



In order to obtain the gain of the other V-band, 2Ai 
A, we must replace A 2 i— > \\. 



The peak gain of the S-band of Eq. (fSJ) , for Ai + A 2 = A, 
is expressed by 



flf 



46 A 2 (Al - A?) 



AiA 2 



{ [£>! (c lC ° - c 2 c°) + ^b] [D 2 (c lC ° - c 2 c°) + E 2 b] } 5 (B2) 



with £>i = 2z/(2^ - A), £> 2 = 2z/(2z; + A), ^ = (c?^ - (Ai - v) 2 ) (44 - \\ + 3v 2 - 2\ 2 v) - Ab 2 (c§) , £ 2 = 

( c?c _ (Aa + ^2) (c?c _ A 2 + 3^ + 2AiJ/) _ 4&2 (c?) 2_ 

Instead, for A 2 — Ai = A, the expression in Eq. ( |B2| ) is valid provided the new coefficients, denoted 
by a prime, are used instead of the unprimed ones, D[ = —2v(2v — A), D 2 = —2v(2v + A), E[ = 

- (44 - (X 1 + v) 2 ) (44 - A 2 + 3v 2 - 2\ 2 v) + Ab 2 (4f, E 2 = (44 - (A 2 + vf) (44 - A 2 + 3v 2 - 2Aii/) + 

46 2 (c?) 2 . 
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